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The single-freeze-out model of hadron production is used to describe the particle 
ratios and the transverse-momentum spectra from RHIC. The emphasis is put on the 
new measurements done at the highest beam energy of ^/sNN—'^^^ GeV. An overall 
very good agreement is found between the data and the model predictions. The data 
for different centrality windows are analyzed separately. A simple scaling of the two 
expansion parameters of the model with the centrality is found. Interestingly, this 
scaling turns out to be equivalent to the scaling of hadron production with the 
number of wounded nucleons. 
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1 Introduction 



Ultra-relativistic heavy-ion collisions offer a unique possibility to produce hot 
and dense hadronic matter in laboratory conditions. With the new results 
presently coming from RHIC [1-3] we have entered an exciting era when var- 
ious model predictions and theoretical ideas may be confronted directly with 
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the large amount of diverse and accurate data. One of the key issues in the 
field is the question of the degree of thermalization of the hadronic matter 
produced in heavy-ion collisions. From the outset of the heavy-ion physics 
the thermal models have been useful and successful to a large degree in the 
description of data [4-17] (for a recent review see [18]). 

In this paper we analyze the hadron production at mid-rapidity measured at 
RHIC at the full beam energy of y^s^N— 200 GeV within the single- freeze- out 
model of Ref. [19 21]. The new element of this work is the application of the 
approach to data at various values of the centrality parameter, c (or the impact 
parameter, 6), which allows us for a closer look at the geometry and flow of the 
hadronic system formed in non-central collisions. The obtained systematics for 
the geometric parameters as functions of c is very interesting, as it conforms 
to the wounded- nucleon scaling [22-24]. We find that the invariant time at 
freezeout, r, and the transverse size, Pmax, behave in a very systematic way 
for not-too-large values of c. In particular pmax ^ {^~^)^ r/p^g,^ c::^ const. 

The single-freeze-out model was formulated in Ref. [19-21] where it was suc- 
cessfully used to describe the transverse-momentum spectra of hadrons mea- 
sured at the then-available lower beam energy of ^/sJ^^130 GeV. The model 
combines a typical thermal approach, used to study the ratios of hadron abun- 
dances [4-11] with a hydrodynamic expansion including in a natural way the 
longitudinal and transverse flow. A characteristic and important feature of the 
model is the complete treatment of the resonances in both the calculation of 
the ratios of hadron yields and the analysis of the spectra. The model and 
the data at ^/sNN=^30 GeV, including strange-quark particles such as 0, A, 
S, as well as X*, have been found to be in a surprising agreement [20]. 
A comparison of the model predictions with the data for particle abundances 
and transverse-momentum spectra, collected at the maximum RHIC energy 
of 200 GeV and presented in detail in this paper, bring further evidence for 
rapid thermalization, or a statistical nature, of the hadronic matter produced 
at mid-rapidity in heavy-ion collisions. Moreover, with our simple approach 
we form a basis for the explanation of most of the soft features of the hadron 
production observed at RHIC. 

The paper is organized as follows: In the next Section we outline the main 
assumptions of the model, bringing up for completeness some of the basic 
material of Refs. [19-21]. In Sect. 3 we study the ratios of hadron yields 
at ^/sNN = 200 GeV and find the thermodynamic parameters characteriz- 
ing the freeze-out. In Sect. 4 we calculate the transverse-momentum spectra. 
The calculations are done for different centrality windows, such as defined by 
the experiment. Naturally, since the geometry of the hot system formed in 
the collision depends on the centrality, we find a dependence of the geomet- 
ric/expansion parameters on c. Interestingly, this dependence refiects quite 
accurately, for not too large c, the scaling of hadron production with the num- 
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ber of wounded nucleons. 



2 The single-freeze-out model 



In our previous work [12,19-21] we have shown that with the complete treat- 
ment of hadronic resonances, the distinction between the traditionally con- 
sidered two freeze-outs, the chemical and the thermal (or kinetic) one (see 
the following discussion), is not necessary. At least for the RHIC energies one 
can achieve a very good explanation of the soft part of the hadronic data by 
assuming a single freeze-out which takes place at the universal temperature 

Tchem ^Icin = T. (1) 

Clearly, a theoretical description of the whole duration of the collision process 
of the two ultra-relativistic heavy ions is difficult, since different degrees of 
freedom are important at various stages of the evolution. Consequently, one is 
forced to describe each stage within a different theoretical framework. In our 
approach we use only hadronic degrees of freedom and thus concentrate on the 
latest stages of the evolution of the system, namely when the hadrons, at some 
instant, cease to interact and stream freely to the detectors. This moment is 
called the freeze-out. Admittedly, this is a far-reaching simplification. Gener- 
ally speaking, the freeze-out in itself may be a complicated process, involving 
duration in time, and a hierarchy where different kinds of particles and dif- 
ferent reactions switch off at different times. In general, by kinetic arguments, 
one expects that reactions with lower cross sections switch off at higher den- 
sities/temperatures, while those with larger cross sections last longer. Since, 
in most cases, the elastic cross sections are larger than the inelastic ones (a 
counter example, though, is the pp interaction), one expects that the inelastic 
(or chemical) freeze-out occurs earlier than the thermal freeze-out. On these 
grounds one may also argue that strange (or charmed [25]) particles decouple 
earlier than other hadrons. With such a picture in mind one may "separate" 
the freeze-out process into a series of more specific freeze-outs associated to 
particular reaction channels. In fact, the works of Heinz and collaborators [26] 
introduced the concept of the chemical and the thermal (or kinetic) freeze-outs, 
with Tchem > Tkin- At the chemical freeze-out the inelastic interactions between 
the hadrons cease, and the chemical composition of the system is fixed. Later, 
only the elastic interactions are effective, leading to further cooling of the sys- 
tem. At the thermal freeze-out even the elastic interactions become ineffective 
and the hadrons are completely decoupled from each other. 

One should bare in mind, however, that the scales in the system responsible 
for the thermalization and expansion are of similar order. To be more precise. 
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for expanding systems the criterion for the thermal freeze-out to occur is [27] 



(2) 



with 



exp 



(3) 



With our parameterization for the expansion profile we obtain Texp — 10 GeV . 
while Ref. [27] gives r^oii ^ lOGeV"^ at T ^ 165 MeV. Thus the discussion 
on the validity of our approach (early thermal freeze-out) as opposed to ap- 
proaches having later thermal freeze-out becomes a discussion concerning the 
precise value of the factor in the criterion of Eq. (2). It is possible that al- 
ready at the moment of hadronization the system is so dilute and possesses 
such a strong flow that not only inelastic, but even the elastic processes can 
be neglected. Thus our main assumption of Eq. (1) may be, in the least, a 
good starting point. At any case, assumptions of a phenomenological model 
like ours are verified a posteriori by the experiment. As we have shown in 
Refs. [19-21] and will also demonstrate in this work for the full RHIC energy 
of a/snn = 200 GeV, we are in a rather comfortable situation. The single- 
freeze-out hypothesis, when combined with expansion, offers a very economic 
description of the data, with very few parameters (two universal thermal pa- 
rameters and two geometric parameters for each centrality bin). 

An experimental argument in favor of a short time between the chemical and 
kinetic freeze-outs has been provided by the observation of resonances, such as 
the iC* or p, by the STAR Collaboration [28-30], in the invariant-mass spectra. 
In the thermal approach the abundance of these resonances is compatible with 
a high freeze-out temperature, around 160 MeV. On the other hand, if the 
elastic scattering processes {K — n for the or tt — tt for the p) took place 
for a much longer time, such that the system would cool off significantly, the 
abundance of the resonances would be much lower. We should point out that 
the scenario of a single freeze-out, put forward in Ref. [19], is compatible with 
the sudden hadronization scenario of Ref. [31]. 

Thus, our starting point is the hypothesis (1). Similarly to other thermal 
models of hadron production, we fix the value of T, and also the value of the 
baryon chemical potential /i^, by fitting the ratios of hadron abundances. The 
results of such an analysis will be presented in the next Section. 

The assumption of the single freezeout allows us to calculate uniformly, i.e. 
within the same model, both the ratios of hadron abundances and the trans- 
verse-momentum spectra. In order to do this, however, we need to make a 
choice of the freeze-out hypersurface {i.e. a three-dimensional volume in the 
four-dimensional space-time) and of the four- velocity field of expansion at the 
freeze-out. A priori, since we do not consider dynamics of the earlier stages 
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here and only make an educated guess, many choices are possible and for 
each particular case one may check the validity of the proposed form of the 
hypersurface and flow by confronting the output of the model calculations 
to the data. Our choice has been made in the spirit of the hydrodynamic 
calculations of Ref. [32-41], and is defined by the condition 



r = 2 _ ^2 _ ^2 _ ^2 ^ const. (4) 
To make the transverse size of the fire-cylinder. 



P = V^^ + i' (5) 
finite, we impose the condition 

P < Pmax. (6) 

In addition, we assume that the four- velocity of the hydrodynamic expansion 
at freeze-out is proportional to the coordinate (Hubble-like expansion), 

T T \ t t t J 

Such a form of the flow at freeze-out, as well as the fact that t and co- 
ordinates are not limited and appear in the boost-invariant combination in 
Eq. (4), imply that our model is boost-invariant. In view of the recent data 
delivered by BRAHMS [42], this assumption turns out to be justifled for the 
description of particle production in the rapidity range — 1 < y < 1, where 
the variation of the observed particle multiplicities is moderate. 



We note that the freeze-out hypersurface and the flow are jointly controlled 
by the parameters r and Pmax- Hence, our model has altogether only four 
parameters (dependent on the colliding energy and the centrality parameter) , 
namely, the two thermodynamic parameters, T and /i^, and the two expansion 
(geometric) parameters, r and Pmax- The thermodynamic parameters turn 
out to depend on the centrality bin very weakly, reflecting the very weak 
observed dependence of the particle ratios. On the other hand, the geometric 
parameters do depend on c and their dependence, obtained from independent 
fitting, assumes a very intuitive form (c/. Sect. 6). 

Other choices for the freezeout hypersurface and expansion have been inves- 
tigated in the literature [43], with a recent study made in Ref. [44], where a 
class of parameterizations has been tested with a limited number of resonances 
included. 

All resonances from the Particle Data Tables [45] have been included in our 
study. The method and necessary formulas that describe the decays in cas- 
cades can be found in Ref. [21]. We stress the important role of including the 
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high-lying states, which although thermally suppressed, are increasingly nu- 
merous according to the Hagedorn hypothesis [46-49]. Moreover, the inclusion 
of resonances is crucial for the description of the p^-spectra, as it increases the 
logarithmic slope as if lowering the temperature [12]. This observation works 
in favor of the hypothesis (1). 





Model 


Experiment 


Fitted thermal parameters 


T [MeV] 


165.6±4.5 




liB [MeV] 


28.5±3.7 




lis [MeV] 


6.9 




Ai/ [MeV] 


-0.9 






0.2 




Ratios used in the thermal analysis 


TT /tt'^ 


1.009 ±0.003 


1.025 ± 0.006 ± 0.018 [50] (0 - 12%) 
1.02 ± 0.02 ± 0.10 [51] (0 - 5%) 


K-/K+ 


0.939 ± 0.008 


0.95 ± 0.03 ± 0.03 [50] (0 - 12%) 
0.92 ± 0.03 ± 0.10 [51] (0 - 5%) 


P/P 


0.74 ± 0.04 


0.73 ± 0.02 ± 0.03 [50] (0 - 12%) 
0.70 ± 0.04 ± 0.10 [51] (0 - 5%) 
0.78 ± 0.05 [52] (0 - 5%) 




0.104 ±0.010 


0.083 ± 0.015 [53] (0 - 5%) 


K'/tt- 


0.174 ±0.001 


0.156 ±0.020 [53] (0 - 5%) 


n/h- X 10^ 


0.990 ±0.120 


0.887 ±0.111 ±0.133 [54] (0- 10%) 


n/h- X 10^ 


0.900 ±0.124 


0.935 ± 0.105 ± 0.140 [54] (0 - 10%) 



Table 1 

The model fit to the particle ratios measured at RHIC at ^/sNN — 200 GeV. Only 
stable hadrons (with respect to the strong interactions) are included. The experi- 
mental and the theoretical pion yields are corrected for the weak decays [28,53]. The 
values in the brackets denote the centrality class. The errors in the model predic- 
tions reflect the full errors of the experimental data points used in the method. 
Note the very good quality of the fit, reflected by the low value of P^r degree of 
freedom. 
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3 Particle ratios 



For boost-invariant systems the ratios of hadron abundances at mid-rapidity 
are equal to the ratios of hadron densities (for a more detailed discussion of 
this point, which includes the effects of the resonance decays, see Ref. [21]). Us- 
ing this property, we can determine thermodynamic parameters of our model 
directly from the study of the ratios of hadron abundances. We stress that 
we have only two independent thermodynamic parameters in our model: the 
temperature, T, and the baryon chemical potential, /llb- The other two chem- 
ical potentials (/i^, /i/) are fixed by the conditions that the initial strangeness 
of the system is zero and that the ratio of the baryon number to the electric 
charge is the same as in the colliding nuclei. This assumption could be lifted 
at the expense of carrying an extra parameter. However, in practice this is 
irrelevant, since the value of /x/ turns out to be very small and the effects of 
the isospin violation can be neglected. 



In Table 1 we show the results of our fit to the RHIC data collected at 's/snn= 
200 GeV. We include here only the stable hadrons (with respect to strong in- 
teractions) and use the most central data. The values of T and fiB shown in 
Table 1 have been fitted with the method to the experimental data listed 
in the third column. We stress a very good quality of the fit, refiected by the 
low value of P^r degree of freedom. An interesting feature of our fit is the 
stable value of the temperature, when compared to other fits done at lower 
energies. The optimum value of T is practically the same as that found in 
the analysis of the RHIC data at -s/snn^ 130 GeV, and (within errors) it 
is the same as the temperature found at the top-energy collisions at CERN 
SPS [55]. ^ On the other hand we observe a clear drop of the baryon chemical 
potential, from 41 MeV found at y^s^N— 130 GeV down to 29 MeV found in 
the present calculation. A smaller value of the baryon chemical potential is a 
simple consequence of the observed increase of the p/p ratio with the beam en- 
ergy. At mid-rapidity, the ratio asymptotes to unity with the increasing beam 
energy. We note the values of parameters shown in Table 1 are compatible 
within the statistical errors wit the results of Braun-Munzinger, Redlich, and 
Stachel [18], who obtain Tchem = 177 ± 7 MeV and /i^ = 29 ± 6 MeV. 



The similarities of the SPS and RHIC spectra are discussed in Ref. [56]. 
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P X [GeV] 

Fig. 1. The transverse- momentum spectra of pions, protons, and kaons measured by 
BRAHMS (diamonds), PHENIX (triangles), PHOBOS (squares), and STAR (stars) 
for the most central events in the Au + Au collisions at ^/snn — 200 GeV. The 
preliminary data are used in the form as compiled in Ref. [59]. The curves represent 
the result of the fit obtained in the single freeze-out model. Full feeding from weak 
decays is included in the calculation. 

4 Spectra 



Let us turn now to the study of the transverse-momentum spectra of hadrons. 
They are calculated with the help of the Cooper-Frye [57,58] formula, 

-p^=/p-dE„/(p^„), (8) 



8 



where the distribution function / includes the products of the resonance decays 
(more details on the construction of / in the presence of the resonance decays 
have been given in Refs. [20,21]). The element of the hypersurface, dS^, is 
defined as 

dx^ dx~^ 

dT.^ = e^a(3^ — da\\ da^dcf), (9) 

where x^ — x^ — r^^ x'^ — Ty^ x^ — and e^c^^^ is the Levi-Civita tensor. 
Introducing a convenient parameterization [37]: 

t = r cosh ay cosh = r sinh ay cosh 

Tx = T sinh a±_ cos 0, '^y — ^ sinh a±_ sin 0, (10) 

one finds 

(iS^(x) = u^{x) sinh(a^)cosh(a^) da^da\\d(j), (11) 

hence the four-vectors rfS^ and are parallel. This feature is special for our 
modeling of the freeze-out, see Eqs. (4) and (7), and allows us to represent the 
spectra in a compact form given by Eq. (8). With the use of the parameteri- 
zation (10) we further rewrite Eq. (8) in the form 

+00 Pmax/r 27r 

— — — ^ J ^^11 J sinha^rf (sinha^) J d^p - u f {p - u) ^ (12) 

^ -oo 

where 

p ' u — mxCOshQf||COshax — Pi. cos ^ sinha^- (13) 

We note that a\i is the rapidity of the fluid element, Vz = r^/t = tanha||, and 
a_i describes the transverse size of the system, p — rsmha_i. 

The maximum and average transverse-flow parameter is given in our model 
by the equations 

/o^"^^^ pdp . ^ 



Since the ratios of hadron abundances depend very weakly on the centrality of 
the collision, we regard our thermodynamic parameters T and /i^ as universal. 
On the other hand, the two geometric parameters may be different for each 
centrality class. From the obvious reasons we expect that r and pmax are larger 
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Fig. 2. Our fit to the transverse-momentum spectra of both the negative and positive 
hadrons measured by BRAHMS [42] for most central events for the Au+Au cohisions 
at y/sNN — 200 GeV. Fuh feeding from weak decays is included in the model 
calculation. 



for more central events or, in other words, they are decreasing functions of the 
centrality, c, defined by the formula 



c = 



C^tot 



(15) 



where b is the impact parameter and atot is the total inelastic cross section 
for the colliding gold nuclei (for a discussion of the validity of this simple 
geometric relation see Ref. [65]). For various centrality bins we fit the geometric 



10 



parameters with the help of the method, i.e. we minimize the expression 

X (^,Pmax) = 2^ ^ 5 (16) 

n=l 

where D^"^ is the n-th measured value of the transverse-momentum spectrum, 
D^°^^^(t, pmax) is the corresponding value calculated in our model, and cr^ is 
the error. Observing the differences between the points measured by different 
experimental groups we conclude that the systematic errors are much larger 
than the statistical errors. Therefore, we assume that is determined by the 
systematic uncertainties and take cr^ = D^^'p/IO.^ 



5 Most central events 



In Fig 1 we show the result of our fit to the transverse-momentum spectra of 
TT", p and for the most central events measured by the four RHIC experi- 
ments ^ (exceptionally, the data of PHOBOS are for the mixtures of negative 
and positive hadrons). In this case, the centrality class of the BRAHMS and 
PHOBOS data is 10%, whereas the centrality class of the PHENIX and STAR 
data is 5%. Our model curves for tt", p and K~ ^ shown in Fig. 1, are simulta- 
neous fits to the four experimental spectra available for each hadron species. 
The fitted values of the geometric parameters are in this case r = 7.58 fm and 
Pmax = 7.27 fm (see Table 2 for a complete list of the values of the geomet- 
ric parameters found for different experimental situations). We observe that 
these numbers are similar to those found in our previous study of the RHIC 
data collected at a lower energy of s/snn = 130 GeV, r = 7.66 fm and pmax 
= 6.69 fm [19,20]. This weak dependence of our geometric parameters on the 
energy of the colliding ions is reminiscent of the very weak energy dependence 
of the measured transverse HBT radii i?side and i?out [61-64]. 

In Fig. 2 we show our fit to the most central BRAHMS data only [42]. In 
this case the fitted values of the geometric parameters are r = 7.68 fm and 
Pmax = 7.46 fm. These values are consistent with the previous values obtained 

^We note that the method was also used to determine the values of the two 
thermodynamic parameters in Sect. 3. Since the thermodynamic parameters are 
sensitive only to the ratios of hadron abundances our fitting of the four parameters 
r, /i^, T, and Pmax can be done in two steps: first we determine T and /i^, and 
having fixed these two values we determine later r and Pmax- One could follow also 
a different path by fitting simultaneously the four parameters to the experimental 
spectra. 

^ Most of the data shown in our figures were digitized from the available plots with 
preliminary experimental results. 
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Fig. 3. Our fit to the transverse-momentum spectra of antiprotons measured by the 
four RHIC experimental groups for the most central events of the Au + Au collisions 
at ^/sNN — 200 GeV. The preliminary data are used in the form as compiled in Ref. 
[59]. The solid curve describes the model result with full feeding of antiprotons from 
the weak decays, whereas the dashed line describes the model calculation without 
any contributions from the weak decays. 

from the "mixed" input. In Figs. 1 and 2 we can see that our model curves 
describe the data very well in the large range of p± stretching up to about 
2 GeV. 

A closer inspection of Figs. 1 and 2 shows that our model curve tends to 
overestimate the spectra of antiprotons at small values of p±. To study the 
possible origin of this effect, in Fig. 3 we show in detail the antiproton spectra 
measured by different experimental groups. One can see that the experimental 
spectra differ substantially from each other and, in practice, one cannot fit 
all four groups of the experimental points simultaneously. The most likely 
reason for such differences are difficulties connected with the correct estimate 
of the role of the weak decays, or with fixing the normalization. Our model 
calculations shown in Fig. 3 indicate that the effect of the weak decays on the 
antiproton spectrum is indeed important. The solid line represents our result 
obtained in the case where the antiprotons are fed by the weak decays. This 
curve may be treated as the upper limit for the spectrum. On the other hand, 
the dashed line represents our result with no feeding included, and may be 
treated as the lower limit. We can see that most of the experimental points. 
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Fig. 4. Our fit to the transverse-momentum spectra of both the negative and positive 
hadrons measured by PHENIX [60] for the Au + Au coUisions at ^/sNN = 200 GeV. 
The eight sets of points and the fitted curves correspond to the eight centrahty bins 
given in Table 2. The two thermal parameters have uniform values for all cases, 
while the geometric parameters are adjusted independently for each centrality bin, 
with results displayed in Table 2. Full feeding from weak decays is included in the 
model calculation. 



except for those of PHOBOS, lie in the band enclosed with our two curves. 
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T»T> A TTTV /TO 

BRAHMS 


10 


7.68 d= 0.19 


7.46 =b 0.05 


0.52 zb 0.01 


0.70 zb 0.01 


STAR 





- 5 


9.74 ± 1.57 


7.74 zb 0.68 


0.45 zb 0.08 


0.62 zb 0.09 




5 - 


- 10 


8.69 ± 1.39 


7.18 zb 0.64 


0.47 zb 0.08 


0.64 zb 0.09 




10 


- 20 


8.12 zb 1.31 


6.44 zb 0.57 


0.45 zb 0.08 


0.62 zb 0.10 




20 


- 30 


7.24 zb 1.18 


5.57 zb 0.50 


0.44 ± 0.08 


0.61 zb 0.10 




30 


-40 


7.07 =b 1.17 


4.63 zb 0.39 


0.39 zb 0.08 


0.55 ± 0.10 




40 


- 50 


6.38 =b 1.02 


3.91 zb 0.33 


0.37 zb 0.07 


0.52 zb 0.09 




50 


- 60 


6.19 ± 1.09 


3.25 ± 0.28 


0.32 ± 0.07 


0.46 ± 0.10 




70 


- 80 


r" AC) 1 o "1 

5.48 lb 0.81 


A f\f^ 1 "1 

4.03 zb 0.10 


0.43 zb 0.06 


0.59 zb 0.07 


PHENIX 





- 5 


7.86 lb 0.38 


7.15 zb 0.13 


0.50 zb 0.02 


0.67 zb 0.02 




20 


- 30 


6.14 zb 0.32 


5.62 zb 0.11 


0.50 zb 0.02 


0.68 lb 0.03 




30 


- 40 


5.73 =b 0.16 


4.95 ± 0.05 


0.48 ± 0.01 


0.65 zb 0.01 




40 


-50 


4.75 ± 0.28 


3.96 ± 0.09 


0.47 ± 0.03 


0.64 zb 0.03 




50 


-60 


3.91 =b 0.23 


3.12 zb 0.07 


0.45 zb 0.03 


0.62 zb 0.03 




60 


-70 


3.67 =b 0.12 


2.67 zb 0.03 


0.42 zb 0.01 


0.59 zb 0.02 




70 


-80 


3.09 lb 0.11 


2.02 zb 0.02 


0.39 zb 0.01 


0.55 zb 0.02 




80 


-91 


2.76 ± 0.20 


1/10 1 r\ 

1.43 ± 0.03 


n oo 1 r\ r\o 

0.32 ± 0.03 


0.46 ± 0.03 


PHENIX + STAR 





-5 


7.96 ± 0.41 


7.15 ± 0.17 


0.50 ± 0.03 


0.67 ± 0.03 




20 


-30 


6.08 ± 0.32 


5.53 ±0.13 


0.50 ± 0.03 


0.67 ± 0.03 




30 


-40 


5.59 ± 0.30 


4.75 ±0.11 


0.48 ± 0.03 


0.65 ± 0.03 




40 


-50 


4.95 ±0.27 


4.03 ±0.10 


0.46 ± 0.03 


0.63 ± 0.03 



Table 2 



The values of the optimum geometric parameters, r and Pmax, and the corresponding 
average and maximum transverse flow velocities, {j3^) and (3^^^ of Eq. (14), obtained 
from our analysis of the data collected by different experimental groups at different 
values of the centrality. All results are for the beam energy of ^/sNN = 200 GeV. 
The model calculation for STAR includes the corrections for the weak decays in a 
way used by the STAR Collaboration, i.e., the feeding of the pions from the decays 
of A is excluded. The model calculations for BRAHMS and PHENIX include fuU 
feeding from the weak decays. 
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A n 



□ K" STAR 200 GeV ' ^ PHENIX200GeV 




0. 5 1 1. 5 2 2. 5 3 0. 5 1 1. 5 2 2. 5 3 

Pi [GeV] 

Fig. 5. Our joint fit to the transverse-momentum spectra from PHENIX [60] 
and STAR [53]: pions (triangles), kaons (squares), and antiprotons measured by 
STAR (open symbols) and PHENIX (filled symbols) for the Au + Au collisions at 
y^SNN — 200 GeV. Full feeding from weak decays is included in the model calcula- 
tion. 

6 Centrality dependence of the geometric parameters 



Having at our disposal the data collected in different centrality windows we 
may determine the dependence of the expansion parameters on c, or, equiv- 
alently, on the impact parameter b. In Fig. 4 we show the PHENIX data 
together with our model curves. Both the data and the model spectra are 
plotted for eight different centrality classes. The most central collision corre- 
spond to the centrality class c = 0-5%, whereas the most peripheral collisions 
included in Fig. 4 are for the centrality class c = 80-91%. The values of the 
fitted expansion parameters for all eight cases are given in Table 2. 

The geometric parameters obtained from the fit to the STAR data turn out to 
be noticeably larger compared to the parameters obtained from the analysis of 
the PHENIX data. In order to see the origin of this effect, in Fig. 5 we superim- 
pose the spectra from STAR and PHENIX for four centrality classes common 
for the two experiments. We can see that the spectra measured by STAR (open 
symbols) are lying slightly above the spectra measured by PHENIX (filled 
symbols). As a consequence of this behavior, our normalization-controlling 
parameter, r, turns out to be larger for STAR. We have performed a simulta- 
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neous fit to the PHENIX and STAR data. The model curves for this case are 
also shown in Fig. 5 (solid, dashed, and long-dashed curves, respectively for 
pions, kaons, and antiprotons) . The fitted values of the expansion parameters 
T and pmax are listed in the lowest entry of Table 2. 




0.2 0.4 0.6 0.8 1 

c 



Fig. 6. Centrality dependence of the geometric parameters, r and Pmax, as extracted 
from the analysis of the PHENIX data only. 




0.2 0.4 ^ 0.6 0.8 1 



Fig. 7. Centrality dependence of the geometric parameters r and Pmax as extracted 
from the analysis of the STAR data only. 

An interesting behavior can be found when the values of r and p^ax are plotted 
versus c (see Figs. 6 and 7). The parameter p^ax exhibits almost ideal linear 
dependence on c in the full range < c < 1. Similarly, the parameter r 
decreases linearly with growing c in the range < c < 0.5. For larger values of 
c the drop of r is slightly weaker, although the data of PHENIX show a rather 
strong tendency to line up. The approximate linear scaling of both r and p^ax 
with c has an important physical significance. The hadron multiplicities in our 
model are functions of and Pmax/"^, see Eq. (12), hence the linear dependence 
of T and pmax on c means that the ratio Pmax/"^ is independent of c and we 
are left with the dependence on only. Consequently, using Eq. (15) we find 
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Fig. 8. The number of wounded nucleons, w{b) (solid line) and the approximating 
function w{0){l—c{b))^ (dashed line), plotted as functions of the impact parameter b. 
Since the multiplicity of hadrons produced in our model is proportional to (1 — c(6))^ 
at moderate values of c, the model conforms to the wounded-nucleon scaling. 

that the total multiplicity in the central region (mid-rapidity) obtained in our 
model should, to a good approximation, scales in the following way 



Nmodelib) = N^odelib = O) (l - c{b)f = N^odel{b = 0) 1 . (17) 



This type of the behavior may be compared to the 6-dependence of hadron 
production following from the wounded nucleon model [22] . For the symmetric 
case, i.e. when two identical nuclei collide, the number of the wounded nucleons 
is given by the formula 



where A is the mass number, is the inelastic nucleon-nucleon cross section, 
and Ta is the nucleon-nucleus thickness function. Similarly, the total inelastic 
cross section is given by the formula 



where Taa (b) is the nucleus-nucleus thickness function. Using the values 
A = 197 and crin=40 mb, and assuming the Woods-Saxon distribution for 




w 




(18) 




(19) 
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the nuclear densities (with the standard choice of the parameters) we find 

w{b) ^ w{0){l - c{b))\ (20) 

Since the wounded-nucleon model assumes the scaling of hadron production 
with the number of the wounded nucleons, the total multiplicity of hadrons 
considered in this model scales also as (1 — c(6))^, hence it is compatible with 
the scaling obtained in our approach, cf. Figs. 6 and 7. The quality of the 
simple approximation (20) is shown in Fig. 8. This result means that our 
model, with the geometric parameters fit to reproduce the data at various 
centralities, conforms with high accuracy to the wounded-nucleon scaling for 
not-too-large values of c (say, c < 0.5). 

We end this Section with a technical remark. In the analysis of this Section 
we have ignored any possible effects of the azimuthal deformation. Clearly, 
the emitting source formed in a non-central collision will not be azimuthally- 
symmetric, which, by the way, is seen in the HBT measurements [66]. However, 
it turns out that the effects of this deformation are tiny for the p^-spectra 
averaged over the azimuthal angle [67], thus can be dropped in the analysis 
of the present paper. 



0. 1 i 




0.5 1 1.5 2 2.5 

- m [GeV] 



Fig. 9. Our fit to the transverse-momentum spectra of ft and Cl, compared measured 
by the four RHIC experimental groups for the most central events of the Au + Au 
collisions at ^/snn = 200 GeV. The preliminary data are taken from Ref. [54]. The 
values of the geometric parameters used in the model calculation are r = 8.3 fm 
and pmax = 7.1 fm. 
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7 Spectrum of the Q baryons 



Finally, we wish to present the p± spectra for the Q and Q. Similarly to 
the case of the lower energy of ^/s^n = 130 GeV [20], our model is in full 
agreement with the preliminary data of Ref. [54], which can be seen in Fig. 9. 
In the model calculation we have used the thermal parameters from Table 1 
and the geometric parameters r = 8.3 fm and p^ax = 7.1 fm, which roughly 
corresponds to centrality — 10% in the STAR experiment (c/. Table 2). 



8 Conclusions 

It is clear from the results presented in this paper that the single freeze- 
out model works well for the particle ratios and the transverse-momentum 
spectra for the full RHIC energy of y^sJvw = 200 GeV. This confirms our 
results published previously for ^/sNN = 130 GeV [19-21]. The model provides 
an economic framework to parameterize the data, and supports the use of 
the thermal approach to heavy-ion physics. In addition, we have shown that 
the model works equally well for the central, as well as for the non-central 
collisions. The transverse-size parameter, Pmax, drops linearly practically for 
all values of c, while the ratio of pmax to r remains almost constant for not- 
too-large values of c. The systematics of the dependence of the geometric/flow 
parameters on the centrality is most interesting. For moderate values of the 
parameter c, say c < 0.5, the analysis of the non-central collisions displays 
scaling in accordance to the wounded-nucleon model. 
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